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The stationary distribution of passive tracers chaotically 
advected by a two-dimensional large-scale flow is investigated. 
The tracer field is forced by resetting the value of the tracer 
in certain localised regions. This problem is mathematically 
equivalent to advection in open flows and results in a fractal 
tracer structure. The spectral exponent of the tracer field is 
different from that for a passive tracer with the usual additive 
forcing (the so called Batchelor spectrum) and is related to 
the fractal dimension of the set of points that have never 
visited the forcing regions. We illustrate this behaviour by 
considering a time-periodic flow whose effect is equivalent to 
a simple two-dimensional area-preserving map. We also show 
that similar structure in the tracer field is found when the 
flow is aperiodic in time. 



I. INTRODUCTION 

It is well known that time-dependent flows with simple 
spatial structure produce complicated filamental patterns 
in advected passive scalars through the phenomenon of 
chaotic advection jjj. Such patterns have been demon- 
strated in laboratory experiments Q and have also 
been simulated numerically e.g. || (in the context of 
laboratory experiments) and |]|7j] (in the context of the 
distribution of long-lived chemical species in the atmo- 
sphere). Similar structures can be observed in satellite 
images of oceanic plankton or sea surface temperature 
||. Although the two or three dimensional distribution 
of atmospheric chemical species is not accessible to high- 
resolution observations, aircraft measurements along one- 
dimensional transects give evidence of filaments and ir- 
regular structure with strong spatial gradients |)fj|fl|l0f . 
The natural way to characterise these complex spatial 
distributions is by investigating scaling properties of sta- 
tistical quantities such as Fourier power spectra or struc- 
ture functions. 

Using concepts from the field of chaotic dynamical sys- 
tems has proved to be useful in understanding the mech- 
anism of mixing of passive scalars (hereafter we shall use 
the term 'tracer' to describe an advected passive scalar) 
in fluid flows with a smooth velocity field and has led to 
the concept of chaotic advection 0] (Lagrangian chaos). 
It has been recognised that the motion of fluid elements is 
typically chaotic even in very simple non-turbulent flows. 
Thus mixing appears as a result of the sensitivity of the 
particle trajectories to their initial position. By stretch- 
ing and folding fluid elements chaotic advection generates 



finer and finer structures and tracer variance cascades 
down to small scales until it is dissipated by molecular 
diffusion and the system ends up in a perfectly mixed 
homogeneous state. 

In order to maintain a nontrivial stationary state it 
is necessary to apply some kind of forcing in the tracer 
evolution equation. As we will show in this paper the 
specific form of this forcing has an important effect on the 
statistical properties of the resulting tracer distribution. 

The usual way to maintain tracer contrast (in numer- 
ical experiments and theory) is by a large-scale additive 
forcing representing sources and sinks of the tracer. For 
this case, theory predicts a power spectrum T(k) ~ fc _1 
for the tracer - the so called Batchelor spectrum - in 
the viscous-convective range, that is between the small- 
est characteristic length scale of the velocity field (L) and 
the characteristic length scale of diffusion (Id)- Vulpiani 
has shown Jl^] that the Batchelor scaling is not restricted 
to the viscous-convective range of turbulent flows but it 
is also valid for non-turbulent (e.g. time-periodic) flows 
if the Lagrangian motion of the fluid elements is chaotic. 
The spectrum has been verified in recent numerical 
simulations 1 13 15. However, attempts to realise a fc _1 
spectrum in laboratory experiments have been inconclu- 
sivegQ. 

In a numerical study |lq |, Pierrehumbert considered a 
different type of forcing - a resetting forcing - in which, as 
fluid parcels entered certain localised regions of the flow, 
tracer concentrations were reset to values associated with 
the relevant region. This was motivated by giving exam- 
ples of a number of experimental settings where a reset- 
ting forcing seems to be more realistic than an additive 
one. Such examples are: i.) thermal convection - where 
the temperature of the fluid parcels is set to the temper- 
ature of the hot (lower) or cold (upper) boundary when 
they cross the thin thermal-boundary layers; ii.) advec- 
tion of dye concentrations - where dye is introduced by 
diffusion from a solid source and maintained at the sat- 
uration concentration in a diffusive boundary layer. We 
note also that the 'resetting' forcing could be relevant 
for chemical tracers (e.g. ozone or nitrous oxide) in the 
lower stratosphere, where polar and equatorial regions 
could play the role of the forcing regions. 

In pj| it was found that the power spectrum of the 
tracer was different from the expected with spectral 
exponent between —1 and —2. However, no precise the- 
oretical explanation for the deviation from the Batchelor 
spectrum was given. 
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Tracer spectra steeper than /c _1 were observed in re- 
cent experiments by Williams et al. fT^] . The injection of 
tracer (dye) in these experiments also seems to resemble 
the 'resetting' forcing. 

In this article, we will show that the deviation from 



the Batchelor spectrum reported in 18 is a direct con- 
sequence of the 'resetting' forcing. Using a dynamicals 
system approach we show that this type of forcing leads 
to the appearance of fractal patterns in the tracer distri- 
bution (similar to the ones characteristic of advection in 
open flows 



that implies a spectral slope steeper 
than —1 (and dependent on the fractal dimension of these 
patterns). In numerical simulations we model advection 
by a simple two-dimensional area-preserving map repre- 
senting a time-periodic velocity field. The effect of ape- 
riodicity of the flow is also investigated, by considering a 
stochastic version of the map. 



II. THE MODEL 

The tracer concentration field c(r, t) in an incompress- 
ible flow is governed by the equation: 



+ v(r, t) ■ Vc = T + DAc 

at 



(1) 



where v(r,i) is the velocity field of the flow. We as- 
sume that v(r, t) is smooth in space with only large scale 
structure (non-turbulent) and time-dependent. Further- 
more, the velocity field is assumed to be independent 
of the tracer concentration (passive tracers). For sim- 
plicity, we will consider two-dimensional flows but we 
expect that similar considerations apply to large-scale 
flows in three dimensions. Two-dimensional flows with 
large scale velocity structure are also relevant for the 
stratosphere, where the transport is mainly along two- 
dimensional isentropic surfaces [2l|p2] |. 

The term T on the right-hand side of (Q) represents 
the forcing. In case of 'additive' forcing T = S(r), rep- 
resenting sources (S > 0) and sinks (S < 0) with a large 
scale structure. The 'resetting' forcing considered in this 
paper can be represented mathematically as a limiting 
case of a relaxation to a concentration profile cq (r) 



T = a(r)(c - c (r)) 
where the rate of relaxation a(r) is 



a(r) 



oo if r G T>i 
otherwise 



(2) 



(3) 



(The case with constant relaxation rate (a(r) = ao) rep- 
resenting some kind of homogeneous chemical or biolog- 
ical activity, was investigated in p3| , p4[ .) Furthermore 
we assume, for simplicity, that Co(rj is constant in each 
forcing domain, cq(t) = C\ for reDj. Thus the concen- 
tration in each fluid element is set to C, as it enters the 



resetting domain T>i. We will study the properties of the 
concentration field outside the forcing domains. 

We note, that tracer advection in a closed flow with the 
resetting forcing applied is mathematically equivalent to 
tracer advection in open flows with time-dependence re- 
stricted to a finite mixing region, ejj. in the wake of a 
cylinder [ ^9|pC| ]. For example, in |2(J] tracer concentra- 
tions were imposed upstream to take one value in part of 
the inflow and a second, different value on the remainder. 
The effect of advection through the wake of the cylinder 
gave a complicated pattern of tracer concentration down- 
stream. In our case, the resetting domains play the role 
of the upstream condition in the open flow problem, since 
fluid parcels emerge from those domains with their value 
of tracer concentration specified. (The resetting domains 
also mimic the downstream escape in the open flow prob- 
lem, since when fluid parcels enter those domains their 
previous value of tracer concentration is immediately for- 
gotten due to the infinitely fast relaxation.) Although in 
p0| the authors noted that the boundary between differ- 
ent values of tracer concentration had a fractal structure, 
the implication for the tracer spectrum was not consid- 
ered. 

The second term on the right hand side of (|]J) repre- 
sents diffusion, which is assumed to be weak. To make 
a more precise statement let us denote by r the charac- 
teristic time of the decay in time of the number of fluid 
particles that have not visited the forcing regions. We as- 
sume that the length scale of forcing (L) and the diffusive 
scale (Id ~ (-D/A) 1 / 2 ) are well separated so that 



1 , L 

T < - log — 

A Id 
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where A is the rate at which chaotic advection generates 
small scale structures, i.e. the absolute value of the most 
negative Lyapunov exponent of the chaotic advection. 
(In two-dimensional incompressible flows this is equal to 
the positive Lyapunov exponent, that is the rate of the 
exponential separation of nearby fluid elements Q.) The 
right hand side of (ji|) is the typical timescale needed for 
the large scale tracer variance to be dissipated by dif- 
fusion [^2| that is much larger then the one needed for 
resetting. Thus, in this case, the fluctuations in tracer 
structure is primarily affected by the resetting and not 
by diffusion. Therefore we neglect diffusion in the follow- 
ings (D = 0) but expect to obtain the same structure as 
it would be for nonzero diffusivity on scales larger then 
Id- 

In the absence of diffusion the concentration is con- 
served along the trajectories of fluid elements (except for 
the resetting) and it is convenient to study the problem 
([!]) in a Lagrangian picture. The concentration for a fluid 
particle depends only on which resetting domain has last 
been visited by that particle. Thus, in order to construct 
the asymptotic tracer distribution we follow the back- 
ward trajectories by integrating 



dr 



(•5) 
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backwards in time from each point until the trajectory 
enters one of the resetting domains and assign the concen- 
tration corresponding to that domain to the point from 
where the trajectory started. The concentration distribu- 
tion along a line segment can be obtained in the same way 
without having to calculate the whole two-dimensional 
field. This allows us to reach high resolution easily for 
one-dimensional transects. 

Advection by a two-dimensional incompressible time- 
periodic velocity field over one time period is equivalent 
to the action of a two-dimensional area-preserving map. 
In general, the explicit form of this map cannot be ob- 
tained analytically from the velocity field. For numerical 
convenience, in the detailed analysis of the tracer field 
that follows we represent advection by a map M. given 
by, 



x n+ i — x n + a sin (2ny n + 4>) mod 1, 
y n+ i = y n + asm (2irx„ + 4>) mod 1, 



(6) 



where x n and y n denotes the coordinates of a fluid par- 
ticle at discrete times t = nT, n — 0, 1, ... and T is the 
period of the flow. It is known that the orbits obtained 
by such iterations are typically chaotic. We note, that 
the action of the map (o) is equivalent to advection for 
one period by a piecewise steady flow given by periodic 
alternation of two sinusoidal shear flows in the x and y 
direction, respectively. 

One characteristic of two-dimensional area-preserving 
maps is the existence of quasiperiodic orbits that form 
KAM tori acting as transport barriers |2(| separating 
different chaotic regions of the phase space. In order to 
avoid such barriers we use a value of the parameter a in 
(^J) that produces a single ergodic chaotic region without 
visible KAM tori. An alternative approach (used in |18|) 
is to consider as a random variable 4> n 6 [0, 2-7r] taking 
on different values in each step of the iteration. We will 
first consider the deterministic case (a = 0.6, <j) n =0), 
and will then discuss the effect of the random parameter, 
representing aperiodic flows. 

Since advection is represented by a map the position 
of fluid parcels is followed explictly only at discrete times 
separated by interval T. We therefore follow Picrrchum- 
bert ( jl8]]) in resetting only at those discrete times and 
also in considering resetting domains of the form of two 
parallel stripes of width e. As we shall see later, vary- 
ing e, and hence the size of the domains, has a strong 
effect on the resulting tracer distribution. [Note that 
every specification of the resetting domains in the time- 
continuous problem defined by (gj), with time-periodic 
flow, (|J) and (|J), corresponds some resetting domain in 
the time-discrete problem. The converse, that every re- 
setting domain in the time-discrete problem corresponds 
to some resetting domain in the time-continuous prob- 
lem, applies provided the slight generalization is allowed 
that a is a function of time as well space in the time- 
continuous problem.] 



III. RESULTS AND DISCUSSION 

To illustrate the characteristic features of the tracer 
field forced by 'resetting' we first consider a time-periodic 
velocity field that is a kinematic model of a meandering 
jet flow psfl . The resetting domains are taken to be par- 
allel to the mean flow below and above the centre of the 
jet, respectively. In Fig.l a snapshot of the concentration 
field (obtained from the backward iteration described in 
the previous section) is shown for this meandering jet 
flow. The periodicity of the flow implies a periodic tracer 
field for long times. In the absence of barriers outside the 
forcing regions almost all trajectories visit the forcing re- 
gions in a finite time. Thus, the concentration field has 
only two possible values (Ci - 'black' and Ci - 'white'). 
The tracer field is formed by elongated filaments with 
a Cantor set like fractal structure transverse to the fila- 
ments (see Fig. lb). The width of the filaments varies in a 
broad range of scales showing that the forcing maintains 
an imperfectly mixed stationary state. Similar structure 
can be observed in the tracer field generated by the map 
(I) (Fig. 2). 

The boundary between the two regions of different con- 
centration value is in fact the boundary between the 
basins of attraction of the two resetting domains for 
the time-reversed advection dynamics. The structure of 
boundaries between basins of attraction of different at- 
tractors is a well known problem in the field of chaotic 
dynamical systems p6|-p8|. Complicated fractal struc- 
ture of basin boundaries appears as a consequence of 
non-attracting chaotic sets, usually associated with the 
phenomenon of transient chaos |2{| . 

Let us consider the set S of points corresponding to 
orbits that never visit the forcing domains 



S = {r\M n (r) £ UjD^Vn = -oo,oo}, 



(7) 



If the map M. is the same in each iteration the set S is, 
by definition, invariant under the iterations of the map, 
Ai(S) = S. Due to the assumed ergodicity of the advec- 
tion, S is necessarily a set of measure zero. Such non- 
attracting invariant sets in chaotic systems (chaotic sad- 
dles) typically have a complex fractal geometry p7| , f28|] . 
In conclusion, the set of forcing domains selects a sub- 
set S of the existing chaotic trajectories that defines an 
invariant fractal set. 

Fig. 3 shows the number of backwards iterations n 
needed to reach either of the resetting domains from 
points lying along a line segment. The boundary between 
the regions with concentration C\ and C2 is formed by 
points for which n is infinite and coincides with the part 
of the unstable manifold of S - to be precise the union 
of the parts of the unstable manifolds of each point in S 
that lie between that point and the first resetting domain 
encountered along the unstable manifold. The results 
of estimating the fractal dimension D of the boundary 
along linear transects using the so called box counting 
algorithm are shown in Fig. 4. 
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For a given flow, the geometry of the resetting regions 
- by selecting an invariant chaotic set - determine the 
fractal dimension D. The dependence of D on the size 
of the resetting stripes e is shown in Fig. 5. As the size 
of the resetting regions increases the fractal dimension 
of the set corresponding to orbits avoiding these regions 
decreases. If the resetting domains occupy the whole 
system, the distribution is formed by a finite number of 
discontinuities (D = 0). (In fact, the fractal dimension 
D may reach zero even if the forcing domains do not 
cover the whole chaotic region.) In the limit of very small 
resetting domains (e — > 0) the set of orbits avoiding them 
increases and S becomes space filling {D —> 1). 

The fractal dimension D of the unstable manifold of 
the chaotic saddle can be related to dynamical properties 
of the system by a variant of the so called Kaplan- Yorke 
formula |3(]] corresponding to non-attracting sets |n| 

D=l-j, (8) 

where A is the Lyapunov exponent on the chaotic saddle 
and k is the escape rate, defined by 

A{nT) = A e- KnT (9) 

where A(nT) is the area occupied by points that have not 
visited the resetting regions by time t = nT, (k = 1/r). 
In our context k might more appropriately be called the 
'resetting rate'. (More precisely, (||) gives the so called 
information dimension pl"[ D% of the invariant measure 
on the saddle, which can be somewhat smaller then the 
box counting dimension). If the non- uniformity of the 
stretching is not too strong the Lyapunov exponent cor- 
responding to any subset of trajectories is close to the 
Lyapunov exponent of the full advection dynamics. In 
our case this means that A has only a weak dependence 
on e and the fractal dimension depends on the structure 
of the resetting regions through the 'resetting rate' k. 

In general, there is no simple formula for the reset- 
ting rate since it depends on the detailed geometry of 
the resetting domains and the structure of the chaotic 
trajectories. In the time-discrete version of the problem, 
however, one can have a situation where mixing is strong 
enough to redistribute the non-reset particles almost uni- 
formly between two resettings, if T 3> 1/A. Then the 
resetting rate can be approximated by 

« = -ln(l-^A 4 ), (10) 

where Aj is the area occupied by the forcing domain % 
normalised by the area of the whole chaotic region. This 
approximation seems to apply to our numerical example 
(A = 2.75/T), as is shown in Fig. 5, where the fractal di- 
mension of the boundary predicted by (^, |Io| ) is compared 
with the measurements. 

The standard way to characterise such tracer distribu- 
tions is by their Fourier power spectra. In our case the 
distribution along a one-dimensional transect is made up 



by an infinite number of jumps between the two possi- 
ble values, with the jumps located on a fractal set (see 
Fig. 2b). The power spectra of such structure has been 
discussed in |32f| , and it was found that there is a sim- 
ple relationship between the spectral slope (7) and the 
fractal dimension (D) of the set of the discontinuities 

7 = 2-L>. (11) 

A more detailed analysis of the relation between fractal 
dimension and spectral slope in p3| and in pij (done in 
the context of advection by compressible flows modelled 
by dissipative maps with fractal attractors) has shown 
that D in ( |TT| ) is the correlation dimension D2 of the 
chaotic set. Assuming further that multifractality (3^] is 
weak (D2 ~ D\) and using (0) the spectral slope can be 
approximated as 

7«l + p (12) 

where the second term represents the deviation from the 
Batchelor scaling and is proportional with the rate of 
'resetting' k. 

For numerical convenience, instead of direct measure- 
ment of the power spectra we consider the second order 
structure function | p6[ defined by 

S(r) = ((C(r + r) - C(r )) 2 } ~ {8rf (13) 

where (..) denotes averaging over different values of the 
coordinate tq. The scaling exponent £ is known to be 
related to the spectral slope by the simple relationship 
7 =1 + C, if 1< 7 <3 @. 

Structure functions (calculated for different transects 
and their average) are shown in Fig. 6. Fluctuations of 
the local slope at large scales decay fast at smaller scales 
where structure functions calculated from different tran- 
sects tend to be parallel. This can be explained by the 
fact that the statistics of the small scale structures is 
much better (as there are many along a transect) than 
that of the few large scale features. In this sense, spatial 
averaging improves the scaling at large or intermediate 
scales. 

The scaling exponent £ depends on the size of the forc- 
ing domains, taking values between £ = 1, (7 = 2) and 
C = 0, (7 = 1). Fig. 7 shows the value of ( as a function 
of the width, e, of the forcing domains. In order to verify 
the relationship ( |TT| ) in our model the measured value 
of the spectral slope (based on the structure function) is 
compared with the one predicted from the fractal dimen- 
sion measurement, for different values of e (Fig. 7). The 
relationship is confirmed with a very good accuracy. 

We also performed numerical simulations with the 
stochastic version of the map (^) representing aperiodic 
flows. The two-dimensional tracer field (Fig. 8) shows a 
pattern qualitatively similar to the one shown in Fig. 2 
but with a somewhat less regular character. An obvi- 
ous difference is that in this case the tracer field changes 
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with time, although its qualitative appearance remains 
the same. 

In the case of aperiodic flows invariant chaotic sets can- 
not exist. Nevertheless, there may exist a set of chaotic 
orbits that never visit the resetting domains. The set of 
points corresponding to these orbits now wanders from 
one iteration to the next, but exhibits the same type of 
fractal structure as seen in the non-random case for the 
set S. (The resetting time, when calculated for points ly- 
ing along one-dimensional transect, also shows a similar 
structure to that in Fig. 3.) 

Structure functions for the aperiodic case are shown in 
Fig. 9. As a consequence of aperiodicity, the spatially av- 
eraged structure function, measured at fixed time, is now 
time-dependent. Furthermore, its local slope fluctuates 
even at relatively small scales. (Compare with Fig. 6b.) 
Nevertheless, the time-average of these structure func- 
tions exhibits clear scaling (Fig. 9b). (Note, that the slope 
of the temporally averaged structure function is different 
from that of the non-random map.) Similar fluctuations 
of the spectrum have also been observed in the case of 
attractors of dissipative random maps p3] and in exper- 
iments with floating particles advected on the surface of 
an aperiodic flow |3§| ] . In both cases, temporal averaging 
led to a good power law scaling. 

The lack of good scaling for the spectra taken at fixed 
times can be explained by the fact that scaling needs 
a good sampling of the possible stretching histories at 
each scale. Structures contributing to the spectra at a 
certain scale I have been advected for a finite time ti ~ 
(1/A) log(L/Z). When the flow is aperiodic (i.e the map 
is different in each time step) the time ti might not be 
enough for a good sampling of the possible velocity fields 
if I is not small enough. This can be compensated by 
temporal averaging, but obviously cannot be replaced by 
spatial averaging in this case. We conclude, that for the 
random map (or aperiodic flows) the power law scaling 
could be restricted to rather small scales but it is clearly 
present in the temporally averaged sense. This also holds 
for the box-counting fractal dimension. There is a well 
defined fractal dimension in a temporally averaged sense 
that relates to the exponent of the structure function as 
in the case of the non-random map. 

IV. SUMMARY AND CONCLUSIONS 

We have shown that when resetting forcing is applied 
to a passive tracer in a chaotic advection flow, the bound- 
aries between different tracer values have a fractal struc- 
ture. Such structure is closely analogous to that previ- 
ously identified in open flows. As a consequence, when 
'resetting' forcing is applied (or, equivalently, when open 
flows are considered) the spectra is steeper than the fc _1 
predicted by Batchelor for the case of additive tracer forc- 
ing. The deviation can be expressed in terms of charac- 
teristic exponents of the advection and forcing dynam- 



ics. The physical mechanism for the steeper spectra is 
that most of the tracer variance is destroyed by the re- 
setting forcing (or escape, in open flows) before reaching 
the dissipative scale. This effect is clearly visible in the 
tracer distribution (Fig. 2) which contains large smooth 
regions, while strong gradients (and hence diffusive dissi- 
pation rates) are concentrated, in the limit of diffusivity 
tending to zero, on a fractal set of measure zero 

The implication of our results is that the form of the 
tracer forcing can strongly affect the spectral slope of 
the advected tracer. Indeed, our results could explain 
why some previous numerical experiments and most of 
the previous laboratory experiments have failed to repro- 
duce the spectrum. Our numerical example shows 
that this is definitely the case for the numerical results 
presented in [l8| . The tracer forcing used in the experi- 
ment |l7|] also seems to be similar to our considerations. 
Although in (l7j the authors try to explain the unex- 
pected deviation from the Batchelor spectra as an inter- 
mittency effect, they recognise that the deviation is much 
stronger then what one would expect for intermittency 
corrections. This was confirmed recently in |l5fl , where 
the effect of intermittency (non-uniform stretching) on 
the spectra was investigated. Non-space-filling, fractal, 
scalar interface in the Batchelor regime was also detected 
in jet flows |39f| where the open flow type behavior, equiv- 
alent to the resetting, should apply. A jet flow was also 
considered in fl6|| reporting spectra steeper then fc _1 in 
the Batchelor regime. The experimental implementation 
of the purely additive forcing seems to be difficult from 
technical point of view and this could be the reason why 
the Batchelor spectra was never unambiguously observed 
in laboratory experiments. 

Of course, the resetting forcing is an extreme limit. 
In reality, the relaxation to a certain tracer value can- 
not be infinitely fast. We expect, however, that similar 
structures would be observed whenever there is a strong 
spatial inhomogeneity of the relaxation rate a(r), hav- 
ing values much smaller than the Lyapunov exponent, A, 
in some mixing dominated regions and being much larger 
then A in others where the relaxational forcing (resetting) 
would have the main effect. 

The 'resetting' model may also be a useful first approx- 
imation for chemical tracers in the lower stratosphere. 
The observed spectra ]40| , 7 w 1.8, fits into the range of 
our prediction (1 < 7 < 2). Furthermore, cliff-like struc- 
tures resembling the ones seen in our simple model are 
characteristic of stratospheric aircraft data. Nonetheless 
those strucutres are not separated by smooth regions. 
The forcing mechanism for stratospheric chemical trac- 
ers is likely to be far more complex than that considered 
here and needs further investigation. 
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FIG. 1. Stationary tracer field (a) (and blow-up (6)) con- 
structed from 400 x 400 particles by calculating backwards 
trajectories until they reach one of the resetting regions. The 
jet model flow and the parameters used are described in [25] . 
The forcing domains are, y > 2.3 - black and y < —2.3 - white. 
Periodic boundary conditions are used in the x coordinate. 



FIG. 2. Stationary tracer distribution (a) (before reset- 
ting) using 400 x 400 particles corresponding to the map ([]) 
using a = 0.6 and <j) — 0- (b), One-dimensional transect for 
x — 0.75. The resetting domains are stripes parallel to the 
y axis centered at x = 0.25-black and x — 0.75 with width 
e = 0.2 



FIG. 3. Number of iterations needed to reach one of the 
forcing domains vs. initial coordinates calculated for 2000 
points on the line segment shown in Fig. 2b. 



G 



FIG. 4. Box counting measurement of the set formed by 
the discontinuities of the tracer field (calculated with high 
resolution using 65536 points) measured along three different 
segments parallel to the x axis y — 0.1; 0.4; 0.7. The slope of 
the power law fit (dashed line) is D = 0.67. 

FIG. 5. Fractal dimension vs. the width of the forcing 
stripes e. Dashed line shows the prediction based on ( JTc| ) 
and (|^), 1 + ln(l — 2e)/A, using the measured value of the 
Lyapunov exponent, A » 2.75. 

FIG. 6. Structure functions calculated for three different 
transects y — 0.1; 0.4; 0.7 (a.) and structure function averaged 
over 10 equidistant transects parallel to the x axis for e = 0.2 
(6.). The power low fit (dashed line) has a slope £ = 0.325. 

FIG. 7. The spectral slope 7, vs. the width of the forcing 
domains, e, based on the measurement of the structure func- 
tion scaling exponent, f, (diamonds), and values predicted 
from the fractal dimension of the boundary (crosses) . The an- 
alytical prediction, 1 — In (1 — 2e)/A, (A = 2.75) is also shown 
(dashed line). 

FIG. 8. Tracer distribution (before the resetting) calcu- 
lated at 400 x 400 points using to the stochastic version of the 
map ([|) with a = 0.6 and <f> n G [0, 2tt]. The resetting domains 
are stripes parallel to the y axis centered at x = 0.25-black 
and x = 0.75 with width e = 0.2 

FIG. 9. Structure functions for the random case, (a.) 
Spatially averaged structure functions at 5 different itera- 
tions, each calculated by averaging aver 10 different transects, 
(e = 0.2). (b.) Temporal average over spatially averaged struc- 
ture functions calculated from 10 consecutive iterations. The 
slope of the power law fit is C = 0.445. 
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